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We present a study of the behavior of two different figures of merit for quantum correlations, entanglement 
of formation and quantum discord, under quantum channels showing how the former can, counterintuitively, 
be more resilient to such environments spoiling effects. By exploiting strict conservation relations between the 
two measures and imposing necessary constraints on the initial conditions we are able to explicitly show this 
predominance is related to build-up of the system-environment correlations. 
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I. INTRODUCTION 

The understanding of quantum correlations has come a long 
way since Schrodingers first mention of entanglement | J]. It 
is well established the potential advantages that quantum cor- 
related systems offer over classical counterparts. Despite a 
substantial understanding of bipartite entanglement, and more 
generally non-classicality, we still lack a fundamental de- 
scription of quantum correlations. Within a particular figure 
of merit's description there are ordering problems associated 
with mixed states while means to quantify and clas- 

sify multipartite correlations are still in the early stages of be- 
ing developed [4- 6] . Quantum correlations are at the core of 
quantum information protocols, and more recently they have 
been linked explicitly to quantum phase transitions I0] and 
even in how biological systems behave 18Q. As it becomes 
clearer that quantum correlations are central to how physical 
systems behave, a complete understanding of them is rightly 
of great importance. 

When addressing the issues arising when dealing with 
mixed states, one is implicitly asking the question of how 
quantum systems interact with their environments j^UTlll. 
Given the fragility of a quantum state, and the impossibility of 
realizing ideal unitary conditions, understanding the dynam- 
ics of quantum correlations in open systems is pivotal. Until 
the discovery of states that do not violate a Bell inequality, 
and yet are inseparable, non-locality and entanglement where 
thought of as the same resource. Now with the plethora of 
means to assess the quantum nature of a state fl2[ [l3ll we are 
faced with a bigger challenge. 

In more recent years the study has significantly broadened 
from the well established notion of quantum entanglement to 
include a wider description aimed at quantifying the genuine 
non classicality of a system. While theoretical advances have 
largely explored these concepts independently more recent 
studies have asked how they relate to each other in relevant 
physical settings, for example dynamics under local chan- 
nels and propagation along spin chains 1 16]. In many 
of these studies quantum discord Il7l ll8ll (QD) has presented 
itself as the more attractive potential resource. One conclusive 
difference between entanglement and QD shown was the latter 
decays asymptotically under Markovian-type noise for a cer- 
tain class of states, while we would see the onset of entangle- 
ment sudden death [19]. Recently interesting works exploring 



the interplay between entanglement and discord Oil [2(1 H]]] 
have similarly shown discord serving as the more optimal re- 
source. This leads to a natural conclusion that the more gen- 
eral picture for quantumness, as captured by QD, is more re- 
silient to noise than entanglement. We show, however, this is 
not always the case and explore conditions under which the 
different notions of quantumness can be contrasted both qual- 
itatively and quantitatively. 

Our study will explicitly address amplitude and phase 
damping environments, two widely applicable forms of noisy 
quantum processes, although the validity of the results can go 
beyond these examples. Using quantifiers with the same en- 
tropic definitions and ensuring initially the measures coincide 
we are able to explicitly show that entanglement of forma- 
tion i22ll (EoF) can, counterintuitively, be shown to be more 
resilient to the environmental action than QD. By exploiting 
a strict conservation relation between the two measures 112 ill 
we are able to show how this is related to the dynamics of the 
correlations developed between the system and environment. 

The remainder of the paper is organized as follows. In 
Sec. [HI we introduce the tools necessary for a complete under- 
standing of the work. Sec. lIII Al explores what happens under 
the action of phase damping environments. Sec. 1MB I covers 
the case of amplitude damping environments. In Sec.lTVlwe 
address the case of mixed initial states. Finally Sec. [V] we 
summarize our conclusions and give some discussions on the 
results. 



II. PRELIMINARIES 

Open quantum systems have attracted considerable theo- 
retical study and as such there is a significant body of tools 
for us to draw from |9D . The study of bipartite quantum cor- 
relations similarly has amassed a number of quantifiers one 
could exploit depending on the specific questions to be ad- 
dressed UHtHI]. As such, we will restrict ourselves to explor- 
ing two different types of damping that embody a wide range 
of physical settings and utilize measures that allow for reason- 
able comparison on the grounds of operational construction. 
We will first consider phase damping environments. Phase 
damping describes the loss of coherence while keeping the 
energy constant. It can be described by the single qubit Kraus 
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Since the phase damping environment only affects the coher- 
ences, it will provide an informative platform to compara- 
tively study different figures of merit for quantum correla- 
tions. We next consider amplitude damping environments. 
This describes the probability of losing an excitation to the 
environment and thus affects both the coherence and energy 
within a system. It can be described via the single qubit Kraus 
operators B23I1 
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The action of the environments on a two qubit state can then 
be calculated 



Q=Y A (K i ®K j )p(K i ®K j )\ 



(3) 



i,j=0 



One can readily determine an extremal asymmetric damping 
scenario, where only a single qubit is damped while the other 
remains completely unaffected, by simply applying the iden- 
tity instead of the Kraus operators to one of the qubits. 

In order to assess how different figures of merit behave un- 
der such channels great care must be taken both at the level 
of initial states and correlation quantifiers. We consider the 
entanglement of formation (EoF) and quantum discord (QD), 
due to their construction they share the same entropic defi- 
nition. The relationship between EoF and QD has been re- 
cently examined to study the distribution of quantum corre- 
lated states JIElllIlHi]]- Moreover, it is possible to establish 
a strict relationship connecting QD and EoF 112 ill , so that such 
two figures of merit appear to be natural choices for a quanti- 
tative comparison 

We take QD 03 H Q as a measure for general quan- 
tum correlations between any two qubits under study. As 
originally proposed by Ollivier and Zurek, QD can be asso- 
ciated with the difference between two classically equivalent 
versions of mutual information, which measures the total cor- 
relations within a quantum state. For a two-qubit state pab, 
the mutual information is defined as 



I(Pab) = S(Pa) + S(Pb) ~ S(Pab). 



(4) 



Here, S(p) = -Tr[p log 2 p] is the von Neumann entropy of the 
generic state p. Alternatively, one can consider the one-way 
classical correlation [181 



(5) 



J^(PAB) = S(p A )-<Hr £U] (A\B) 9 



where we have introduced < H{ft i }(A\B)=J] i piS(p l A \ B ) as the 
quantum conditional entropy associated with the the post- 
measurement density matrix p l A ^ B = Tr b^iPabM Pi obtained 
upon performing the complete projective measurement {II/} 
on qubit B. QD is thus defined as 



TT =inf [I(pAB)-J^(g A B)] 



(6) 



with the infimum calculated over the set of projectors {ft/} JTvL 
l25ll . Analogously, one can define which is obtained upon 
swapping the roles of A and B. 

On the other hand, our chosen entanglement measure is 
EoF |22|], which quantifies the minimum number of Bell pairs 
needed in order to prepare a copy of the state we are studying. 
For arbitrary two-qubit states, EoF is calculated as 



(7) 



where h(x) = - xlog 2 x - (1 - x)log 2 (l - x) is the binary en- 
tropy function and Q„ is the concurrence of the state Ii22l1. 
The latter, an equally valid entanglement measure, is found 
in terms of the eigenvalues A\ > ^2,3,4 of the matrix pAB(d~y ® 
^pAB^y ® & y) as 



Ci n = max 



0, V^-^ 



i=2 



(8) 



where fr y is the y-Pauli operator. 

Comparing two measures of quantum correlation stands as 
fundamental problem. Even if we consider the same figure of 
merit, e.g. entanglement, we are plagued by ordering prob- 
lems between different measures when dealing with an arbi- 
trary (mixed) state. Such problems become compounded if 
one attempts to draw conclusions between two different fig- 
ures of merit, such as entanglement and discord, for an arbi- 
trarily chosen state. Therefore we will restrict our study to 
initially pure states of the form 



|0> = a|00>+ Vl - or 2 111) 
PF) = or|01>+ Vl -a 2 |10) 



(9) 
(10) 



where we have parameterized the amount of initial correlation 
in the state in terms of the concurrence d„ as 



1+ Ji-cl 



(11) 



As mentioned, when the state at hand is pure, the EoF and 
QD are exactly equivalent, both quantitatively and from an 
operational view point, and are given by the von Neumann 
entropy. This will allow us to have a justifiable starting point 
from which we will be able to compare the measures even 
when the state has become mixed due to the interaction with 
the environments. 



III. DYNAMICS OF QUANTUM CORRELATIONS 
A. Phase Damping 

The study of phase damping provides a fertile area to study 
the behavior of quantum correlations since the nature of the 
noise is one which only affects the quantum coherences and 
will leave the classical correlations intact. Through Eqs. (1) 
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FIG. 1 : EoF black (dashed) and QD red (solid) for a phase damping 
affected pair of initially pure qubits. Each plot is for an increasing 
amount of initial correlation corresponding to a C in = (a) 0.5, and (b) 
1.0. 



and (3) we can readily find the damped density matrix and ex- 
plore the how the different construction of correlation quanti- 
fiers behave under the same type of noisy channel. A further 
remark is in order: it is immaterial which form of initial state 
Eq. © or (TTTTb we consider. This is a consequence of the fact 
the channel only effects the off diagonal terms appearing in the 
density matrix and has no affect on the populations. Hence it 
does not matter which subspace of the Hilbert space we are 
confined to. 

By examining Fig. [T] we find that regardless of the initial 
correlations and of the strength of the channel, QD decreases 
at a faster rate than EoF. This would appear counterintuitive 
considering the broad definition of the the two figures of merit, 
QD by nature includes states that are entangled, and indeed 
one cannot have an entangled state with zero QD. Thus the 
natural answer to such a counterintuitive result is to conclude 
that the two measures are fundamentally incomparable. How- 
ever, given the carefully chosen initial states and figures of 
merit, we see that the two measures are identical in their def- 
inition initially. By taking these conditions into account one 
must reach the conclusion that quantitatively EoF performs 
better than QD under phase damping. 

In order to explain why this is the case we scale the prob- 
lem back so we can make use of the results of Ref. [21]. If we 
consider an asymmetric damping, such that only one qubit in 
the initially correlated pair is damped and without loss of gen- 
erality assuming Eq. (4) as our initial state the total system- 
environment state can be written 

|^>^ = or |000>+ Vl-a 2 (Vl-^|101)+ VI|111>). 

(12) 

This state is a pure tripartite entangled state. In Ref. W21I1 the 
authors proved a strict conservation relation between the EoF 
and QD for any mixed bipartite states by examining the dis- 
tribution of the correlations in a tripartite pure state. As stated 
inJH 

"Given an arbitrary tripartite pure system, the 
sum of all possible bipartite entanglement shared 
with a particular subsystem, as given by the EoF, 
cannot be increased without increasing, by the 
same amount, the sum of all QD shared with this 
same subsystem." 
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FIG. 2: Quantum discord of (a) p AB , (b) p AEk and (c) Pe a e b - Entan- 
glement is zero everywhere except in p AB . Plots for phase damping 
affecting both qubits. See text for discussion 



Hence, 

The asymmetry usually present when dealing with QD is not 
an issue for us as under this noise QD is symmetrical for our 
initial states. If we now examine what happens for Eq. (IT2l) 
we find that while the damped qubit A develops a non-zero 
QD with the environment E, the pair never share any entan- 
glement, i.e. Sae = V A, a. Hence 

S AB = D AB + D AE , (14) 

therefore QD, in this case, can never be larger than EoF. For 
the case of only one of the qubits damped this explains the 
better resilience of EoF to phase damping over QD. As X —> 1 
we see all quantum correlations are lost and we are left with a 
completely classical total state of qubit-environment. 

When both qubits are damped the system exhibits an anal- 
ogous behavior. However we are unable to exploit any strict 
relations like Eq. dT3b . Also as shown in |2g, l27|] EoF and 
QD do not, in general, exhibit monogamic behavior. In Fig. [2] 
we plot the QD for pab, Pae a , and Pe a e b - Here it is clear as 
the correlation shared between the two qubits is reduced due 
to their interaction with the environments there is a non-zero 
QD developed between all the other bipartitions of the com- 
posite system. However, with the exception of pab, the en- 
tanglement in all of these bipartitions is exactly zero and thus 
the entanglement is quantitatively and qualitatively larger than 
discord. Although for ease of calculation we have assumed the 
damping rate in each environment is the same, such a behav- 
ior evidently holds even when the damping parameters for the 
environments are different. 
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FIG. 3: EoF of (a) p A B, (b) Pae and (c) p B E and one-way QD of (d) p A B, (e) Pae and (f) p 5£ for amplitude damping affecting only qubit A. See 
text for discussion. 



B. Amplitude Damping 

We now consider the amplitude damping channel given by 
Eq. 0. If we first consider the case in which only a single 
qubit is damped we can once again use the results of 112 ill 
to contrast quantitatively EoF and QD although care must be 
taken. The form of the initial state can now have an impact 
on the dynamics of quantum correlations due to the changing 
energy within the system. The states resulting from the action 
of the channel on qubit A only is 



l*AoW = tf|ooo>+ Vi-a 2 (Vrloii)+ Vwino))(i5) 

\fa) ABE = a\010)+ Vl-a 2 (Vrl001)+ VWU00>)(.16) 

Due to the asymmetry of QD now we must be strict in the 
subsystem that we perform the measurements on when calcu- 
lating the QD in order to make use of the conservation relation 
Eq. (IT3T) . If we restrict to only measurement on A then the re- 
duced state dynamics for all pairs will be independent of what 
the initial state was. This is simply because by restricting to 
measurements on A the formula for QD does not take into ac- 
count the differences arising by starting from an initial state in 
the single- or two-excitation subspace. 

The qualitative behavior of the system is very different un- 
der amplitude damping. There is clearly an energy exchange 
happening between the system and environment and therefore 
it is intuitive to expect some build of correlations between the 
system and environment. And, in fact, when the qubit A is 
maximally damped all correlations between A and B are lost. 
However we find the reduced state pbe now shares precisely 
the same amount of initial correlations. Hence the amplitude 
damping admits an effective entanglement swapping action 
whereby the correlations are exchanged from A to the envi- 
ronment. This is in contrast to the phase damping case ad- 



dress previously, since there only coherences were diminished 
and there is no exchange of correlations when the system is 
maximally damped. In Fig. [3] we plot the EoF and one-way 
QD, measured on the damped qubit A, between the various 
reduced states of Eq. (fT5l) . The first major difference we no- 
tice is Sae ± 0. As Eq. dT3T) is pure, and making use of the 
conservation relation Eq. ([T3lh we find regardless of y or a 



(17) 



with equality only holding when y = or 1. Hence by the 
conservation relation the EoF in pab will always be greater 
than the QD. In this case we again see entanglement serves as 
a more robust form of quantum correlations over discord. 

If we consider a symmetrized version of QD, sometimes 
referred to as "two-way" QD 112511 . 



D s = Max[£r,£T], 



(18) 



Eq. (TT3T) no longer holds. Although EoF typically performs 
better than QD there are instances when the initial state has 
Ci n - 1 when QD is the better conserved quantity. How- 
ever, making definitive statements comparing the two mea- 
sures is dangerous. Entanglement by construction makes no 
distinction between the subsystems within the state while dis- 
cord does. The asymmetry of QD poses a significant issue if 
we wish to make quantitative comparisons. We are perfectly 
justified in concluding EoF is a preferred quantity when we 
consider amplitude damping if we measure only on qubit A, 
while the "two-way" QD cannot be used to compare the mea- 
sures since it is no longer fully rooted in the same operational 
definition as EoF due to the maximization in Eq. (H~8l) . 

When both qubits are damped these issues become com- 
pounded significantly, as shown in Fig. 01 When both qubits 
are damped, both EoF and QD are now sensitive to the form of 
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(a) (b) (c) 




FIG. 4: EoF (dashed, black), one-way QD (solid, red) for 2-qubits initially in | x F>-state, Eq. JTOl , and EoF (gray, dot-dashed), one-way QD 
(green, dotted) for 2-qubits initially in |0>-state, Eq. (|9j both affected by amplitude damping channels. All insets show the symmetrized version 
of QD Eq. dTST) for the same states. In all plots C in =0.5. 



the initial state. This means that any firm conclusions on the 
comparative behavior is difficult. Clearly there are regions in 
which one measure appears preferable to another, but with the 
lack of any strict relation there is little meaning to the com- 
parison. 



IV. MIXED INITIAL STATES 

For completeness we now turn our attention to mixed initial 
states. Exploring the behavior for a general two qubit mixed 
state is computationally extremely costly and provides little 
insight to our study. To reduce the calculation complexity we 
examine what happens when the initial state of the two qubits 
is in a Werner state 1 28] 

Q W = n\iff)(iff\ + ^^-i. (19) 

where \if/) is a maximally entangled pure Bell state. As our 
starting state is no longer pure we find EoF and QD do not 
coincide initially. In fact, Eq. (|T9l) is separable for rj < 1/3 
while it has a non-zero discord for the whole range of 7] e 
[0, 1). In order to attempt to remove this discrepancy we shall 
focus on a rescaled quantity for the correlations, i.e. G/Gi„ 
(D/Di n ), where G in (D in ) is the initial value of EoF (QD) for 
the undamped system. The result of this means we are asking 
what portion of correlations, as quantified by a given figure of 
merit, are degraded by the noisy environment. 

Fig. \5\ shows the behavior of QD and EoF of Eq. (H~9b for 
phase damping, panels (a) and (b), and amplitude damping, 
panels (c) and (d). In the case of phase damping if the initial 
correlation is not large enough we find the rescaled QD ap- 
pears more robust to noise than EoF, however the larger the 
initial correlation, and hence the less mixed the initial state, 
we see this behavior reversed. Amplitude damping behaves 
quite differently. Here the rescaled QD is almost always out- 
performing EoF, until the state is nearly pure. This behavior 
is explained in 11 1 611 where the ability to propagate nonclas- 
sical correlations was explicitly linked to the probability of 
exchanging an excitation. 



Clearly, quantitatively comparing EoF and QD here is dan- 
gerous. Despite taking a rescaled quantity in an attempt to 
make the measures comparable, the fact that QD and EoF do 
not coincide initially leads to ambiguous results. Before we 
can make a definitive conclusion a deeper understanding of the 
interplay between entanglement and discord for mixed states 
is needed. 



V. CONCLUSIONS AND DISCUSSIONS 

In this paper we have examined how both quantitatively 
and qualitatively different figures of merit for quantum cor- 
relations behave under noisy processes. In order to make 
such comparisons we considered two correlation measures 
that share the same entropic definitions, namely entanglement 



(a) (b) 




FIG. 5: Rescaled EoF (dashed, black) and rescaled QD (solid, red) 
for the 2-qubit mixed state Eq. [\9l under going (a), (b) phase damp- 
ing and (c), (d) amplitude damping. 77 = 0.5 panels (a) and (c), 
77 = 0.95 panels (b) and (d). 
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of formation and quantum discord. These measures are such 
that when the state is pure, both are equivalent and equal to 
the von Neumann entropy. We considered a broad class of 
states that are pure and parameterized by their concurrence. 
Under these circumstances we see it can be meaningful to 
make both qualitative and quantitative comparisons between 
the two measures of quantum correlation. Interestingly we 
have shown the EoF serves as a preferred resource over QD 
when the noisy process affects only the quantum correlations 
and leaves classical correlations and the energies unaffected. 
By exploiting the strict relation between EoF and QD derived 
in 112 ill we were able to show that the reason for EoFs better 
performance was due to the fact that at no point during the 
dynamical evolution did the damped system ever develop en- 
tanglement with the environment, while there was a non-zero 
QD developed. When the noisy process was one such that 
the exchange of excitations between the system and environ- 
ment was allowed, drawing firm conclusions was more diffi- 
cult. The results highlighted another issue arising from the 
asymmetry of QD. When exploring the effects of amplitude 
damping, due the exchange of excitations with the environ- 
ment the QD was sensitive to the form of the initial state and 
on which subsystem we measured. If we consider only one- 
way QD and only a one qubit of the pair damped, then by the 
same reasoning described for phase damping we can see the 
EoF is a more robust quantifier of quantum correlations. How- 
ever, considering the symmetrized QD we can no longer be so 
conclusive. And when we considered both qubits damped we 
lack any tools to make a meaningful comparison. 

Making comparisons between measures has stood as a diffi- 
cult task in quantum theory. Here we consider some necessary 
constraints in order to quantitatively compare the measures: 
the two measures are operationally equivalent and, must coin- 
cide initially. In Il20ll a strict relationship between negativity 
and the geometric measure of quantum discord was proven 
for all two qubit states, N 2 < Qd- If we consider the exact 
same analysis performed in Sec. IIII Al we find this inequal- 
ity is saturated. However, examining the dynamics of the re- 
duced states of pae we still find a non-zero geometric dis- 
cord, while negativity is exactly zero. This makes compar- 
ing the behavior of the two measures meaningless and high- 
lights some recent concerns associated with the geometric dis- 
cord 1 29]. It may be tempting to rely on other distance based 
measures such as the relative entropy of entanglement (E R ) 
and discord (D R ). However, in this case D R is likely to be 
unfairly favored since the set of classical states is a subset of 



the set of separable states, i.e. a state cannot be classical and 
non-separable, while the converse is easily possible. The dif- 
ficultly in attempting to use these measures is more evident 
if we consider a phase damping affected pair of qubits. Since 
the noise only affects the coherences it is fairly easy to see that 
the closest classical state is exactly the closest separable state 
resulting in E R = D R , however there is still QD developed be- 
tween the environment and system while entanglement is ex- 
actly zero. The same situation encountered when considering 
squared negativity and geometric discord. If we remove the 
condition for the states to be initially pure, and instead con- 
sider the Werner states as our initial states, then once again 
we lose the ability to make firm conclusions. This is because 
the two measures are no longer equal due to the initial mixed- 
ness. 

The results raise a number of important questions, princi- 
pally: when can we compare entanglement and discord? Here 
we made a physically motivated approach and explored the 
measures behavior under noise. The conditions considered al- 
lowed us to show how the conservation relation Eq. (IT3l) can 
be used to make a comparative picture of quantum correla- 
tions. Also, some studies have shown QD to perform better 
than entanglement, e.g. Ill4l4l6ll . something that would ap- 
pear intuitively normal given the definition of the quantifiers, 
while here we show the converse can also be true. Indeed, 
if we view quantum correlations as potential information pro- 
cessing resources then it is important to understand what is 
the advantage of exploiting one form of correlation over the 
other. The sole resilience to noise does not imply advantages 
of one potential resource over the others. While discord has 
been shown to be an attractive candidate for information pro- 
cessing (3(1 H2] our results highlight the need for a deeper 
understanding of the nature of quantum correlations at a fun- 
damental level. We expect our results will serve as an ini- 
tial step to determining a broader framework for the under- 
standing and comparison between different figures of merit 
for quantum correlations. 
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